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I. INTRODUCTION 

The topological susceptibility of the ground state of Yang-Mills theories is an important 
parameter for understanding the breaking of the U(l) axial symmetry. It is defined by the 
vacuum expectation value 

X = j d'x{0\T{q{x)q{0)}\0), (1) 

where 

q{x) = -£^,e^'''"'F^^{x)F^M (2) 

is the topological charge density. Lattice is the ideal tool to determine x (fo^' ^ review cf. 
e.g. Ref. i). 

On the lattice a topological charge density operator with the appropriate classical con- 
tinuum limit can be defined as 

1 

^^^^^ ^ ~ 24 X 327r2 S e^^p^Ti [U^Mpa] , (3) 

fiupa=±l 

where U^^, is the product of link variables around a plaquette. Its classical continuum limit 
is 

g^(a;) — . a' q{x) + 0{a') , (4) 

where a is the lattice spacing. We define the lattice bare topological susceptibility from 
the correlation at zero momentum of two q^{x) operators: 

X X 

X^ is connected to x by a non-trivial relationship: the presence of irrelevant operators of 
higher dimension in q^{x) and in the lattice action induces quantum corrections. Eq. 
must be therefore corrected by including a renormalization function Z{l3) 0: 



q^{x) — ^ a'^Z{/3)q{x) + O(a^) . (6) 

Further contributions arise from contact terms, i.e. from the singular limit a; — >^ of Eq. (^. 
They appear as mixings with the trace of the energy-momentum tensor 

T{x) = ^F^Mx) (7) 



and with the unity operator, which are the only available renormalization-group-invariant 
operators with dimension equal or lower than x- Therefore the relationship between the 
lattice and the continuum topological susceptibility takes the form 
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= a'Z\P)x + a'A{P) (T(x))n.p. + + 0(a6) , (8) 



where N.P. denotes taking the non-perturbative part. T{x) is the correct operator to describe 
the mixing P , since it has no anomalous dimension. In order to be consistent with traditional 
notation, we will rescale T(x) to the gluon condensate 6*2, which is defined by 

where &o = lliV/487r^ is the first coefficient of the /^-function. In Eq. (P) the proportionality 
constant is fixed by requiring 

G2 = ^|^f;,f;,(x)^ + o{g'). (10) 

We will then write Eq. (|]) in the form 

X\P) = a'Z\i5)x + aM(/5)G2 + P(/5)(/) + 0{a'') . (11) 

Z(/3), P{(3), and A{f3) are ultraviolet effects, i.e. they have their origin in the ultraviolet 
cut-off- dependent modes. They can be computed in perturbation theory. In the present 
paper, we will show explicitly that Z{(3), P{P) and A{P) (computed non-perturbatively) are 
well approximated by the first few perturbative terms. This fact is far from trivial, since 
in this case the perturbative series are not even Borel-summable (cf. Ref. Moreover, 
in principle they could be affected by non-perturbative contributions; however, there are 
arguments p| implying that non-perturbative effects should not appear in the first few 
terms of the perturbative expansion. 

In Ref. a "heating" method was proposed to estimate Z{P) and -P(/5) directly from 
numerical simulations, without using perturbation theory. This method has already been 
employed to study, on the lattice, the topological properties of 2-d CP^~^ models |]9HT2 



which are toy models enjoying many similarities with QCD. In the present paper we apply 
the heating method to the A-d SU(2) Yang-Mills theory on the lattice, and we show that it 
can also be used to disentangle the contribution of the mixing with the operator T{x). 

The main idea of the method is to exploit the fact that renormalizations are produced by 
short-ranged quantum fiuctuations, at the scale of the cut-off a, whereas physical effects hke 
gluon condensation or topological properties involve distances of the order of the correlation 
length. When using local updating procedures in Monte Carlo simulations, ffuctuations at 
distance / ~ a are soon thermalized, whereas ffuctuations at the scale of the correlation 
length are critically slowed down when approaching the continuum limit. For a standard 
local algorithm, e.g. Metropolis or heat bath, the number of sweeps needed to thermalize 
the ffuctuations at distance ^ should grow proportionally to with 2; ~ 2. Quantities like 
topological charge, involving changes of global properties of the configurations, are expected 
to be even slower to reach equilibrium, experiencing a more severe form of critical slowing 
down. This is suggested by the fact that in the cooling procedure the topological charge 
survives long after the disappearance of the string tension [|T^. In large- lattice CP^~^ 



models the autocorrelation time of the topological susceptibility grows exponentially with 
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respect to the correlation length [O]. A similar phenomenon has been observed in the 2-d 



U(l) gauge model |]T5 



The heating method consists in constructing on the lattice a smooth configuration car- 
rying a definite topological charge Qq, and heating it by a local updating procedure at 
a given value of [3. Short-ranged fiuctuations, contributing to Z{[3) and P{P), are rapidly 
thermalized, while the initial global topological structure is preserved for a much larger num- 
ber of local updatings. This allows us to obtain estimates of Z{j3) when heating instanton 
configurations and of the mixings when heating fiat configurations. 

In particular when heating a fiat configuration (a configuration with zero fields, i.e. with 
all the link variables equal to the identity), which has zero topological charge, we expect 
the production of instantons at the scale of ^ to take place much later than the thermaliza- 
tion of local quantum fiuctuations. The above-mentioned considerations on critical slowing 
down effects lead to the expectation that, when ^ 3> a, the heating procedure on a flat 
conflguration should show the following intermediate stages before reaching full equilibrium: 

(a) Short-ranged fluctuations at / ~ a contributing to P{f3) get thermalized in a number 
of updating sweeps independent of ^. The signal of should show a plateau giving 

(b) When the number of heating sweeps n increases up to oc fluctuations at / ~ ^ 
start to be thermalized, and gluon condensate and its mixing with sets on. This should 
produce an increase in x^- 

(c) Since a more severe form of critical slowing down is expected to affect the topological 
properties, the x^ signal should show a second plateau in which the topological charge is 
still zero and the whole signal is mixing to the identity operator and to the trace of the 
energy-momentum tensor. This can be checked by cooling back the sample of conflgurations 
and controlling that Q = 0. 

(d) Eventually the modes responsible of the topological structure will be thermalized 
and x^ 'w^ill reach its equilibrium value. 

Actually, in the SU(2) lattice gauge theory the situation is complicated by the fact that 
we are forced to work at small ^ (low jS): the term involving x iii Eq. ( pT]) is exponentially 
suppressed with respect to the mixing with the identity operator; therefore it becomes 
rapidly smaller than the errors and is not detectable at larger ^ (large (3). Using the Wilson 
action, the optimal region where to investigate Eq. (|TT|) is around (3 = 2.5. For these values 
of /3, the correlation length as obtained from the square root of the string tension is ^o- ~ 4, 
while from the lowest glueball mass one obtains ~ 2. Of course the relevant correlation 
length depends on the quantity we are studying. For example, since the gluon condensate is 
closely connected to the glueball propagation, the relevant correlation length should be ^g. 
In this case of small ^, the flrst two regimes (a) and (b) in the heating procedure may become 
not clearly distinguishable. However, we will still have a clear intermediate plateau where 
the topological structure is trivial (stage (c)), allowing us to separate the pure topological 
contribution from the mixing terms in Eq. (|11|)- 

In Section H we present our numerical results. They include data from standard Monte 
Carlo simulations and from the heating procedure. In Section ^ we draw our conclusions. 



II. NUMERICAL RESULTS 
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A. Monte Carlo simulations and perturbative results 



We performed standard Monte Carlo simulations on a 12^ lattice using the Wilson action 
and collecting data for over an extended range of /3. We employed the over-heat-bath 
updating procedure [|16[. We also measured the topological susceptibility Xcooi by using the 
cooling method which consists in measuring the topological susceptibility on an en- 



semble of configurations cooled by locally minimizing the actions (starting from equilibrium 
configurations). The topological content of the cooled configurations is measured by using 



The topological susceptibility measured on cooled configurations is seen to 



gradually reach a long plateau. We estimate Xcooi from the plateau measurements. Data for 
and Xcooi are reported in Table |. 



In Eq. (|TTD , A{(5) and -P(/?) can be calculated in perturbation theory following the 

field theory prescriptions. The first few terms of the series are 



where z\ has been calculated finding z\ = —2.1448 [0], 



MP) 



f3^ ■■■ ' 



where bo = 1.874 x 10"^ iUHl, and 



C3 C4 



P' 



+ 



(12) 



(13) 



(14) 



where cg = 2.648 x 10"^ and C4 = 0.700 x 10"^ [0. 

While the first calculated terms of P{P) fit the data well at large (3 (/5 > 4), more terms 
must be included as j3 decreases. Lacking an analytical calculation for these terms, one 
must fit them from the data. In order to estimate them we, fitted data for f3 > 2.8, where 
the whole non-perturbative signal is smaller than the errors. One more term proved to be 
sufficient for a fit with x^/d.o.f ~ 1. We found 



C5 



3.6(4) X 10" 



(15) 



We extrapolate this result to get an estimate of P{P) at /? ~ 2.5. An overall fit to the 
data of Table | gives a value of C5 consistent with (0), and allows to extract the signal of 
dimension 4 in Eq. (0): 



Ai 



+ 



Ai 



0.3 X 10^ 



(16) 



for f3 ~ 2.5. There is no way within this method to separate the term proportional to x from 
the mixing to G2 without both a direct computation of Z{(3) and A{f3) and an independent 
determination of e.g. by cooling. In fact the term proportional to G2 is indistinguishable 
from contributions 0(l//3^) (or higher) in Z'^{f3). 

Evidence for the existence of the mixing to G2 was obtained by comparing different 



definitions of xl [10|. In what follows we will instead obtain the two contributions indepen- 
dently. 
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B. Heating an instanton configuration 



We start from a configuration Co which is an approximate minimum of the lattice action 
and carries a definite topological charge Qq (typically Qq ~ ±1). We heat it by a local 
updating procedure in order to introduce short-ranged fluctuations, taking care to leave the 
background topological structure unchanged. We construct ensembles C„ of many indepen- 
dent configurations obtained by heating the starting configuration Cq for the same number 
n of updating steps, and average the topological charge over C„ at fixed n. Fluctuations of 
length / ~ a should rapidly thermalize, while the topological structure of the initial config- 
uration is left unchanged for a long time. If, for a given /3, we plot = J2x averaged 
over Cn as a function of n, we should observe first a decrease of the signal, originated by the 
onset of Z{i3) during thermalization of the short-ranged modes, followed by a plateau. The 
average of over plateau configurations should be approximately equal to Z{P) Qq. Since 
we do not expect short-ranged fluctuations to be critically slowed down, the starting point 
of the plateau should be independent of /?. 

In order to check that heating does not change the background topological structure of 
the initial configuration, after a given number Uc of heating sweeps we cool the configurations 
(by locally minimizing the action) and verify that the cooled configurations have topological 
charge equal to Qq. 

We remind that the size of our lattice is 12^. As heating procedure we used the heat-bath 
algorithm, which is efficient in updating short-ranged fluctuations, but is severely affected 
by critical slowing down for larger modes. We constructed on the lattice an instanton 
configuration according to the method described in Ref. |2l[]. On a 12^ lattice we found that 



the optimal size of the instanton is p = 4. We performed also a few cooling steps in order 
to smooth over the configuration at the lattice periodic boundary. After this procedure 
we end up with a smooth configuration with Q^ ~ 0.90 (on the lattice Q^ measured on a 
discrete approximation to an instanton is exactly 1 only for very large instantons and in the 
infinite- volume limit). 

In Fig. |l]we plot Q^{Cn)/Qo, where Q^{Cn) is the lattice topological charge Q^ averaged 
over the ensemble Cn- The data in Fig. ^were taken at /3 = 2.5 and (3 = 3.0. We see clearly a 
plateau starting from n = 4 for both values of (3. The check of the stability of the topological 
structure was performed at ric = 5. According to the above-mentioned considerations, the 
value of Q^/Qq at the plateau gives an estimate of Z{P). We repeated this procedure for 
other values of /?. The behavior of Q^(C„)/Qo is always very similar to the case reported in 
Fig. D. The results are presented in Table |I|. 

A fit of the data to a polynomial 



ZiP) = 1 - ±i|^2 + 9 (17) 



2.1448 Z2 

gives Z2 = 0.48(4) with xV^-o.f. ~ 0.3. 

In Ref. 0] Z2 was estimated by comparing Xcooi with the signal of dimension 4, and it 
was determined to be ~ 1.2. However, it is clear from Eq. (^) that the estimate included 
the mixing with G2. By using the new determination of Z{(3), we can extract A{f3) G2 from 
the old data; approximating A{f3) with the first term 62//?^ is consistent with the data and 
gives 
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^ ^ 0.5(2) X 10« . (18) 



C. Heating a flat configuration 

We now proceed to the analysis of the ensembles C„ of configurations obtained by heating 
the flat configuration (using the heat-bath algorithm), for several values of {3. In Figs. ^ and 
^ we plot the average value of as a function of the number n of heating steps respectively 
for [3 = 2.45 and /3 = 2.5. At = 10 for /5 = 2.45 and = 10, 15 for /5 = 2.5 we check by 
cooling that the topological charge is still zero. 

For both values of (3 we observe long plateaus starting from n ~ 10, which are lower 
than the equilibrium values of (see Table |ID, but also definitely higher than the estimate 
of P{P) obtained in Section [11 A| . We repeated this procedure for other values of /3. The 
behavior of x^(Cn) is always similar to the cases plotted in the figures. Data for the quantity 
measured on the plateaus Xpi are reported in Table 

We believe that Xpi contains the mixings both to the identity operator and to the trace 
of the energy-momentum tensor for the following reasons: 

(i) The plateaus observed are long and after cooling back we do not find any topological 
structure. 

(ii) The plateau values of x^ are systematically higher then the values of -P(/3) obtained 
from Eqs. (|14D and (||). 

(iii) At /? ~ 2.5, the correlation length relevant to the gluon condensate is small (it comes 
from the lowest gluebal mass: ^j, ~ 2), therefore at n ~ the fluctuations contributing to 
the gluon condensate start to be thermalized; for n > 10 they could be already approximately 
thermalized. 

(iv) Data are not inconsistent with a first plateau at P{P), followed by a second one; they 
are also not inconsistent with the expected shift of the starting point of the second plateau, 
when f3 is increased from 2.45 to 2.5. However, since the change in correlation length is 
small, these phenomena can not be detected clearly within our error bars. 

The quantity 

should then measure the physical topological susceptibility. Data for Xh are reported in 
Table |. In order to evaluate XhiP), we inserted into Eq. (^) the parametrization for Z{P) 
given by Eq. (^), using the fitted value of Z2- 

III. CONCLUSIONS 

In the previous section we obtained two independent estimates of the topological suscep- 
tibility: Xcooi, by cooling method, and Xh, given by the relationship (|19|). The comparison 
is satisfactory, although Xcooi seems to be systematically lower. This behavior could be 
explained by the fact that underestimates the topological charge content of the cooled 
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configurations (for tlie lattice size we are working witli), as we found out explicitly when we 
constructed an instanton configuration on the lattice (cf. Section IIB ). 



In the limit P ^ oo, the lattice spacing a is given by the two-loop renormalization 
formula 

1 / c \ 51/121 / o 



- j^^fiP) , fm = [j^^'P j exp (^"YY^ /5J • (20) 

In Fig. Ijwe plot Xh/-^L Xcooi/A|^ versus /3. Scaling is quite good within the errors. By 
fitting to a constant we found: 



|| = 3.5(4) X 10^ , ^ 
According to the arguments given in Section ^1 C|, the difference: 



- 3.5(4) X 10^ , ^ = 2.8(2) x 10^ . (21) 



M{(3) ^ xtm-P{f3) (22) 

should be proportional to the gluon condensate. Assuming that the first non-trivial term of 
the perturbative expansion ( Eq. (p!3| ) ) gives a good approximation of the function A{l3), 
we can estimate G2 by 

G2 - Gm = ^ ■ (23) 

In Fig. ^ we plot Gm/ Although the errors are large, the signal is clear. Fitting the 
data to a constant we obtained 

^ = 0.38(6) X 10« . (24) 
Al 

This result is consistent with Eq. (|18]). This estimate can also be compared with an in- 
dependent determination obtained by studying the plaquette operator [^,^: G2I ^\ = 
0.30(2) X 10^ 0. The comparison is again satisfactory. 

We have shown that Eq. (|TT|) is physically meaningful on the lattice, since it separates 
contributions having different physical origin. Mixings with the unity operator and with 
the trace of the energy-momentum tensor are well defined, and we have estimated them. 
Estimates of x and G2 coming from Eq. ( pJ]) and the heating method are consistent with 
those obtained independently by other methods. 

Finally it is pleasant to notice that our results strongly support the feasibility of the 
determination of G2 from the lattice topological susceptibility in full QCD, i.e. in the presence 
of dynamical quarks, as proposed in Ref. ||23 . 
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FIGURES 



FIG. 1. Determination of the multiplicative renormalization Z{(3). Dashed lines indicate the 
value of Z{P) estimated by averaging data on the plateau. 

FIG. 2. vs. the number of updatings when heating a flat configuration at /3 = 2.45. The 
dashed lines indicate the equilibrium value of and the dot-dashed lines the value of P{(3) (with 
the respective errors). The solid line shows the value of Xp\ estimated by averaging data on the 
plateau. The result obtained by cooling the configurations after ric heating sweeps is indicated 
with the symbol ■. 

FIG. 3. x^ vs. the number of updatings when heating a flat configuration at /? = 2.5. The 
dashed lines indicate the equilibrium value of x^ and the dot-dashed lines the value of P{(3) (with 
the respective errors). The solid line shows the value of Xp\ estimated by averaging data on the 
plateau. The results obtained by cooling the configurations after Uc heating sweeps are reported 
with the symbol ■. 

FIG. 4. Plot of Xh/^i (O), Xcooi/Ai ( X ) and Gm /A-i ( ■ ) vs. /3. Data for Xh and Xcooi 
are slightly displaced for sake of readibility. 
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TABLES 



TABLE I. 

X^5 Xcooij obtained by standard Monte Carlo simulations, and Xh vs. (3. The column 
"stat" reports the statistics of the Monte Carlo simulations. 



(3 Stat lO^x^ lO^Xcooi lO^Xh 



2.45 


80k 


3.14(3) 


7.7(5) 


9.8(1.8) 


2.475 


80k 


2.91(3) 


6.5(4) 


7.6(1.2) 


2.5 


60k 


2.69(4) 


4.4(3) 


5.1(1.1) 


2.525 


80k 


2.56(2) 


3.0(3) 


4.6(1.0) 


2.7 


4k 


1.84(6) 








OK 


1 P,C\( A\ 
i.uvy^ j 






2.9 


4k 


1.36(4) 






3.0 


7k 


1.21(3) 






3.25 


4k 


0.93(3) 






3.5 


4k 


0.69(2) 






3.75 


4k 


0.58(2) 






4.0 


3k 


0.45(2) 






4.5 


3k 


0.33(1) 






5.0 


3k 


0.227(8) 






6.0 


3k 


0.127(5) 







TABLE II. Measure of the multiplicative renormalization of , starting from an instanton of 
size p = 4 on a lattice 12"^ and with Qq = 0.90. The estimate of Z{(5) is taken by averaging the 
data in the range of n reported in the column "plateau" . Since data on the plateau are correlated, 
as error we report the typical error of data in the plateau. 



3 Stat plateau 

2A5 2k ^ 0.20(2) 

2.475 3k 4-7 0.22(2) 

2.5 5k 4-7 0.22(1) 
2.55 3k 5-7 0.23(1) 

2.6 2k 4-7 0.25(2) 
2.8 Ik 5-7 0.32(2) 
3.0 0.6k 4-7 0.33(2) 
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TABLE III. Xpi '^s. (3 starting from a flat configuration. Xp\ is estimated by averaging the data 
in the range of n reported in the column "plateau" . As error we report the typical error of data in 
the plateau. 



p 


stat 


plateau 


10^X^1 


2.45 


5k 


10-20 


2.72(6) 


2.475 


8k 


10-15 


2.56(4) 


2.5 


10k 


11-15 


2.44(4) 


2.525 


10k 


10-15 


2.32(4) 
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